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Abstract. Two causes of non-locality inherent in nucleon-nucleus scattering are considered. They are the 
results of two-nucleon antisymmetry of the projectile with each nucleon in the nucleus and the dynamic 
polarization potential representation of channel coupling. For energies ~ 40 — 300 MeV, a g-folding model 
of the optical potential is used to show the influence of the knock-out process that is a result of the 
two-nucleon antisymmetry. To explore the dynamic polarization potential caused by channel coupling, a 
multichannel algebraic scattering model has been used for low-energy scattering. 



PACS. 24.10.-i Nuclear reactions models and methods 
methods - 24.10.Ht Optical and diffraction models 



24.10.Eq Coupled channel and distorted wave 



1 Introduction 

Non-locality is omnipresent in modern physics having rele- 
vance to a quite diverse array of topics. Examples are the 
heavy-element abundance patterns in white dwarf stars, 
the bulk modulus of silicon, the tensionless limit of string 
theory, non-commutative quantum field theory, image re- 
trieval in scanning probe microscopy, electron crystallog- 
raphy, Bose-Einstein condensates, quantum nanostructures, 
high Reynolds number flow and optical recognition of mu- 
sical scores. 

In nuclear physics, the nucleon-based structure of the 
nucleus and the scattering and reactions of a projectile 
(nucleon or nucleus) with a target nucleus, are compli- 
cated many-body problems. They involve non-local inter- 
actions, and one of the major sources of non-locality arises 
from the effects of the Pauli principle. The many-nucleon 
wave function of the system must be antisymmetric with 
respect to interchange of all pairs of nucleons. This leads to 
the interaction of an individual nucleon with the rest of the 
system (the "optical potential" in the case of scattering) 
having both direct and exchange (knock-out) components. 
The exchange components are intrinsically non-local. Be- 
cause of the short-ranged nature of the strong nucleon- 
nucleon (NN) interaction, the exchange amplitudes have 
significant magnitudes; ones comparable to direct scatter- 
ing amplitudes. 

In this paper, we deal specifically with non-localities 
arising in effective interaction methods of modelling elas- 



tic nucleon-nucleus (NA) scattering. The optical poten- 
tial resulting from such theories is complex (when there 
is flux loss from the incident channel), energy dependent, 
and non-local. These characteristics are due in part to the 
complex nature of the effective NN interaction underly- 
ing the process being modified from the free NN force 
because of the presence of the nuclear medium and from 
effects of the Pauli principle [F. A non-local and complex 
optical potential results also from truncating the Hilbert 
space 0. 

The most common approach to analyze (elastic) scat- 
tering data is the phenomcnological optical model (OM) 
in which, usually, the NA interaction potential is chosen 
to have Woods-Saxon (WS) form with complex strengths, 
and central plus spin-orbit terms. Those optical poten- 
tials are local in form with the energy dependence of the 
parameters assumed to account for non-local effects. One 
exception is the specific energy-independent, non-local po- 
tential studied by Perey and Buck [3]. There have been 
other non-local NA OM potentials proposed, but those 
invariably have parameters to be specified in searches for 
good fits to each specific elastic-scattering data set. Such 
phenomenological approaches, and particularly wave func- 
tions derived therefrom, are not reliable. The actual non- 
localities resulting from more precise treatment of Pauli 
principle effects, such as with the exchange (knock-out) 
amplitudes in scattering, and/or of specific accounting for 
channel-coupling processes, lead to specific energy depen- 
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dent non-locality effects that are comparable with the non- 
exchange amplitudes and/or local form corrections to a lo- 
cal optical potential for the channel-coupling processes. In 
part, the problems are hidden by parameter adjustments 
in the approximate models being able to give asymptotic 
phase shifts which lead to reasonable fits to cross-section 
and other data. Of course, there are a number of ways to 
determine a local equivalent potential given a specific non- 
local one. A review of these forms Chapter 11 in Ref. [Tj. 
The results usually are energy dependent and usually in- 
volve differentials. However, to fit data with these poten- 
tials still requires determination of the parameter set in- 
volved. A significant effect of the violation of the Pauli 
principle inherent may persist [4]. Nonetheless a major 
problem with use of any local equivalent form lies with 
just how large the off-diagonal properties of the non-local 
optical potentials are and/or how large, and with what 
interference, are their contributions to scattering. 

For different energy regimes, when one seeks to min- 
imize phenomenology, different models of NA scattering 
are relevant. At energies greater than the excitation ener- 
gies of giant resonances in the target, a g- folding model p] 
built upon the Kerman-McManus-Thaler (KMT) theory 
of scattering |5j define microscopic optical potentials. That 
approach uses complex, energy- and density-dependent 
NN effective interactions. With optical potentials built 
using the Melbourne force [1] as the effective interaction, 
nucleon-nucleus scattering over the entire mass range and 
for energies ~ 40 to over 200 MeV has been analyzed suc- 
cessfully [TJ . In the regime of energies coinciding with the 
excitation values of giant resonances, one needs consider 
their virtual excitation at least as a second-order exchange 
process [6]. However, at energies from to those for which 
the giant resonances are expected to be dominant door- 
way states, and when the target spectrum is discrete, it 
is more appropriate to use a coupled-channel method to 
analyze scattering data. That is especially so whenever 
the discrete spectrum of the target has low- lying collective 
states. However, most coupled-channel methods of solving 
NA scattering problems are based on a collective-model 
formalism and limited to using local form factors. These 
methods approximate non-localities with local effects by 
taking the parameters, used in finding solution of the cou- 
pled equations in coordinate space, to be energy depen- 
dent. But in so doing the Pauli principle is violated with 
serious effects [7]. That can be overcome by using the, 
recently developed [5], multi-channel algebraic scattering 
(MCAS) scheme; applications of which, even starting with 
interactions built from the simplest, local form, collec- 
tive model descriptions, have been quite successful [4,8, 
[9] . This method of solution of coupled-channel scattering 
problems has a number of salient features, notably: 

1. With this method, it is quite straight-forward to define 
specific radial forms of the NA optical potentials which 
account for the effects of all channels considered [10] . 

2. Pauli exclusion effects can be included even when sim- 
ple (local) collective models are used to define the ini- 
tial potential function matrix V cc t(r). 



3. All sub-threshold (bound) states of the compound sys- 
tem can be defined; 

4. There is a procedure by which the centroids and widths 
of all resonances in the scattering can be determined, 
no matter how narrow any resonance may be. 

5. It can be used with any matrix of initial potentials, 
whether they be local or non-local in form. The lat- 
ter will be the case, due to the Pauli principle, when 
detailed nuclear structure is used to define the initial 
matrix of interaction potentials. 

6. S'-matrices and, in fact, complete off-shell T-matrices 
are evaluated; the latter required for analysis of off- 
shell processes such as bremsstrahlung or nucleon cap- 
ture reaction cross sections. 

In the next section, we present some theoretical con- 
siderations, outlining the methods we have used to ana- 
lyze scattering data and specify the effects of non-locality 
in the problems. For the case of low-energy scattering, 
we also identify the non-locality in the optical potential. 
Then, in Sec. Ill, we present and discuss results. Finally, 
in Sec. IV, conclusions are drawn. 



2 Theory considerations 

An elegant formulation of NA optical potentials, given in 
the book by Feshbach [2], has an operator specification 
of the dynamical polarization potential (DPP). The Fesh- 
bach formalism uses projection operators that divide the 
Hilbert space into the channels of a given scattering prob- 
lem that are to be considered explicitly, V, and all other 
channels, Q. Feshbach determined that the Schrodinger 
equation for scattering took the form 

(E - H vv - H VQ [E - Hqq + if]' 1 H QT "j |^+)) = , 

(1) 

where HpQ = VHQ and similarly for Hqp, Hpp, and 
Hqq. Thus the DPP, formally, is a projection of Q space 
effects to an effective operator in V space, i.e. 



AUpp = HpQ [E — Hqq 



ieP 1 H, 



QV 



(2) 



We consider the case that V projects onto the elastic- 
scattering channel, then 

V = \V ga ) {V gs \ ; Q = 1 - \W g .) (W ga \ ; Q \9 g .) = , (3) 



and take as the Hamiltonian, 
H = H Q + V + H A ■ 



Ha |*j> = ej 



(4) 



\<Pj) are the eigenstates of the Hamiltonian for the tar- 
get nucleus and ej are the eigenvalues. Then as only the 
operator V connects Q and V spaces, the DPP is 



AUpp — VpQ [H — Hqq 



1 Vqv . 



(5) 



For simplicity, temporarily, we ignore antisymmetriza- 
tion between the continuum projectile nucleon and all the 
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A nucleons of the target [2]. Then the (A + 1) particle 
states have the form 



|^) = |x; +) )l^) = |x+(0)) \^(l,---,A)) , 
V W) = K) = |x+(0)) \$ gs {l,---,A)) 
and on taking the target ground-state expectation, 



(6) 



E-H 



($ a s\V\<P g 



- <K 



fix 



)J |x+)=0, (7) 



where 



G 



(+) 
QQ 



[E — Hqq + ie]~ 



(8) 



Thus, a structural form for the NA optical potential is 
identified by 



Uqm(E) 



($ gs \v\$ g s) 



= (* B .\v\$ 



y.i I 



as 

AU(E) 



(9) 



Here AU(E) is the DPP which leads to a coordinate space 
optical potential, 

U M{r,r';E)=V l0C (r)5(r-v') + AU(r,r';E) , (10) 

when a local form (V/ oc ) for the elemental ground-state 
interaction is assumed. 

But allowing for the Pauli principle whereby the emer- 
gent nucleon may not be that incident on the target makes 
even the leading term non-local. The Schrodinger equa- 
tions for the relative motion wave functions thus have the 
form (with x^ + H r ) = X ( r >-^0) 

^ 2 + E ) X (+) (r)= J UoM(r,r';E) X (+) (r')dr> . 

'"_ ' ^ (11) 

Here m is the reduced mass, and the energy scale is taken 
with e gs = 



2.1 Low-energy regime and MCAS 

For nucleons of up to 5 MeV scattering from light mass 
targets particularly, cross sections show relatively few (res- 
onance) states which are often widely spaced. Low ener- 
gies also mean few partial waves are important with the 
scattering (both factors keep the scale of the problem of 
using MCAS manageable). Nonetheless, the applicability 
of MCAS is only limited by the computing power and time 
available. 

The essential input to MCAS is a matrix of poten- 
tials that define a specific NA system. The first oper- 
ational requirement then is to find the optimal expan- 
sion of those potentials in separable form. That expan- 
sion is made in terms of sturmian functions [8, 10J with 
the sturmians being generated from the chosen matrix of 
potentials themselves. Then, with separable interactions, 



the Hilbert-Schmidt expansion of amplitudes gives the T- 
matrix also in separable form. 

The MCAS theory [8] treats coupled-channel scatter- 
ing in momentum space, giving solutions of Lippmann- 
Schwinger (LS) integral equations. For each spin parity 
J n of a given NA system, one must consider a set of r 
scattering channels, each with a label c (1 < c < T) where 
each c identifies a set of quantum numbers (all details have 
been published [8]). Then, in partial wave form, the cou- 
pled LS equations define a multichannel T-matrix, 



T 



(p,q;E) = V c 

open 

E 



cc 



V, 



T J,(p,x)- 



, q; E) dx 



closed 



c" = l 



hi, + x 1 



rj-ij 



(x, q\ E) dx 



(12) 



E) for E > e c and E < e c respectively. e c is the 



where /i = 2m/h 2 . There are two summations as the 
open and closed channel components have been separated 
with the wave numbers being k c 

energy threshold at which channel c is open. Note that 
henceforth the J 71 ' superscript is to be understood. Ex- 
pansion of V cc ' in terms of a finite number of sturmians 
yields 

N 

Vcc'iP, ?) ~ E Xcn(p) Vn 1 Xc'n(q) , (13) 

71 = 1 

where Xcn(p) relate to the chosen sturmian functions in 
momentum space and 77 are the associated eigenvalues [8j 
[TO]. 

The Fourier-Bessel transforms of the form factors are 
Xcn(f)- They are also defined in terms of the sturmians 
($cm(r)) by 

r 

Xcn(r) = J2 V <*'( r )*c'n(r) (14) 

c' = l 

if the initial potentials are local in form, and by 
r poo 

Xc P (0 = E/ V cc ,(r,r')$ c , p (r')dr' (15) 
d=i Ja 

if those potentials are non-local. It is possible to develop 
the optical potential in coordinate space with the MCAS 
theory [8lH0]. We will consider only spin zero targets so 
that the channel indices designate the collected angular 
momenta c = Ulh)jl; J] ■ Then, for each conserved total 
angular momentum J and selected parity, the partial wave 
I associated with a zero spin ground state is unique. In the 
following, the index c = 1 then implies the partial wave 
specific to each J™ considered. 

Assuming a local form for the leading term, MCAS 
theory gives for the DPP, 



r 

E 

c,c'=2 



AU 1 , 1 (r,r';E)= V V lc (r) G { C Q J (r, r'; E) V c , x {r') . 



(16) 
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The Green's functions are solutions of 

r 



Gif=G^6 cc ,+ Y: GPVc^a&l, (17) 



c"=2 



where G c is the free Green's function for each individ- 
ual channel c. Historically, finding these solutions has been 
difficult. However the MCAS development allows the anal- 
ogous definition pTjQj], 



N 



AU 11 (r,r';E) = £ Xln (r) [A(E)] m , Xln ,(r') , (18) 



to 24.0 MeV, and for a set of targets ranging from 27 A1 
to 208 Pb. Good agreement was found between results ob- 
tained with those parameters and data taken with those 
targets. 

Though the PB interaction was not used for scattering 
from 12 C, we adopt it in this study to display its non-local 
features. The PB prescription is energy independent and 
so the shape of that non- locality is constant. The results 
may not be the PB potential that gives a best fit to scat- 
tering data from 12 C. But as we only seek to make quali- 
tative comparison with the MCAS results it suffices to use 
the parameter values initially defined [5] . 



where 



A(E)= V -G^(E) 



G { Q \e) 



open oo 



E 



x 2 + ie 



1 dx 



closed -qo „ 

. J2 f Xcn{x ) Xcn 'S x) #,u- 



\ Jo hi + x 2 



(19) 



2.1.1 The Perey-Buck (PB) non-local potential 

The PB potential [3] has the energy- independent, non- 
local, form 



2.2 The g-folding model for medium energies 



For energies well above the particle-emission threshold and 
so coinciding with a continuum of states in the target, 
multi-nucleon scattering reductions of the formal Fesh- 
bach theory of scattering are relevant. It is useful to fol- 
low the KMT scheme [5] in which it is assumed that only 
pairwise interactions between the projectile and individual 
target nucleons are important. With the projectile tagged 
by the subscript 0, and target nucleons tagged similarly 
by i, the essential interaction can be written 



V 



A 



Av 01 



(23) 



r'|) 



1 



[y/irf3 NL } 



- f r~r' 



' Y-Y' 



where, with s = |(r + r'), 
v(s) = V NL 



I —Rnl \ 

A a NL J 



1 -1 



(20) 



(21) 



As shown in the appendix of Ref . [3] , this (reduced) form 
can be expanded in partial waves with radial multipoles 
being 

1 



—9t(ry) = 3 v(i[r + r'}) 
rr 1 V^Pnl 



2rr' 



(22) 



However, there are parameters involved and so this poten- 
tial is phenomenological. Nonetheless, the parameters of 
this interaction, intended for use at low energies and with 
nuclei of medium and heavy mass, were obtained solely 
by fitting the differential cross sections for 7.0 and 14.5 
MeV neutrons elastically scattered from 208 Pb. Those pa- 
rameter values then were used [5] in analyses of a set of 
neutron scattering observables for a range of energy, 0.4 



With whatever is chosen for the free NN interaction «oi , 
the KMT theory allows for multiple pairwise interactions 
so that JVJV ^-matrices, solutions of NN LS equations, are 
required in applications. However, such an approach does 
not account for the ways in which those interactions are 
influenced by the nuclear medium in which the two nu- 
cleons interact. Experience [T] has shown that the TV TV 
^-matrices, solutions of the Brueckner-Bethe-Goldstone 
(BBG) equations for infinite nuclear matter [13j . can well 
approximate those (many-body) corrections. Using the NN 
(^-matrices in a first order KMT approach involves an as- 
sumption that the important terms in Q-space, for the 
incident energies concerned, are excitations in which a 
particle is promoted to a state in an infinite-matter sys- 
tem with Fermi momentum related to the density of the 
nucleus from where it came. Thus, the effective JVJV in- 
teraction, besides being energy dependent and complex, 
also is density dependent pQ. The optical potentials that 
result on folding these interactions with density matrices 
of the target nucleus, will be both complex and energy- 
dependent. Further, by virtue of the Pauli principle at 
least, they are non-local. 

The g-matrices can be used directly within a momen- 
tum space formulation of the NA optical potential [T3] 
or by forming an effective NA interaction in coordinate 
space [I]. Our approach is with the latter and requires a 
mapping of the actual NN (7-matrices into the coordinate 
space forms usable in the program DWBA98 [15] . 
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In coordinate space, such non-local optical potentials 
can be written formally 

Uom = UoM(rx,r 2 ;E) = Ug M {r x ;E) + [/^(n,r 2 ;£) 

(24) 

where 

USM{ r u E ) = ^ZCn5{r l -r 2 ) / ip* n (s) v D (r ls ) tp n (s) d 3 s 
= S(n-r 2 ) I p{s)v D {r ls )cPs (25) 



which allows expansion of the scattering amplitudes in the 
form of weighted two-nucleon elements since a,j m (l) \&jm) 
in Eq. (|28[) is independent of coordinate '1'. Thus 



involving then an integration over the nuclear density, and 
the exchange term which is totally nonlocal, 



](n<Pn{ r l) v Ex(ri 2 )<p n (r 2 ) . (26) 



vd , vex are combinations of the components of the ef- 
fective NN interactions, (p n {r) are single-nucleon bound 
states, and ( n are bound-state shell occupancies. More 
generally, the latter are one-body density matrix elements 
(OBDME). It is not easy to specify actual radial values of 
these non-local potentials for graphing purposes. The pro- 
gram DWBA98 [15] does not form such in evaluation of 
solutions of the partial wave Schrodinger equations. The 
operator structure and strong density dependence of the 
effective interactions deduced from the NN (/-matrices do 
not lend themselves easily to create multipole expansions 
needed to specify Uom^i, r 2 ) for each partial wave. In 
finding solutions of the integro-diffcrcntial equations, the 
program DWBA98 uses expansions involving particle-hole 
expectation values of the effective interactions PQ. How- 
ever, the effects of the inherent non-locality are clearly ev- 
ident from calculated cross sections and spin observables 
as we show later. They are also important in predictions 
of cross-section and spin-observable data from inelastic 
scattering of nucleons from nuclei. 

Cross sections for inelastic proton scattering have been 
evaluated using a microscopic DWA theory of the pro- 
cesses pp. In that theory, the transition amplitudes for 
nucleon inelastic scattering from a nuclear target have the 
form pp (!?j 4 Af i = !Pj 4 M- i (l---A)) 

T = T M t M iV ' V{Qsc) 

= (x^k^l^M,! Ag cff (0,l) 

xA m {\ X ( +\m) l#W},(27) 

where f2 sc is the scattering angle and Aoi is the two- 
nucleon state antisymmctrization operator. The nuclear 
transition is from a state \ JiMi) to a state \JfMf) and 
the projectile has spin projections v before, and v' after, 
the collision. The incoming and outgoing distorted waves 
are specified by an d they have relative momenta of lq 
and k D respectively. The development proceeds by using 
a cofactor expansion of the target states, 

\*jm) = ^ E \Vim{l)) a jm (l) \*j M ) , (28) 



x (xlHKO) (^ 2TO2 (1)| goff(0,l) 



oi {|xS +) (k,0)) |^ imi (l)>} 



I 



(-l)°^ mi) (jU 2 ™i -m 2 \I -N) 



31,3 
X 



x(x^(k o 0) (<p jama (l)\ g cff (0,l) 

Ai{|xl +) (k,0)) |p iimi (l))},(29) 

where reduction of the structure factor to (transition) OB- 
DME [Sj^j^ij f° r angular momentum transfer values 

/ is given in detail elsewhere pp. 

The effective interactions g c s(0, 1) used in the folding 
to get the optical potentials have also been used as the 
transition operators effecting the excitations. As with the 
generation of the elastic-scattering optical potentials from 
which the distorted waves are generated, antisymmetry 
of the projectile with each individual bound nucleon is 
treated exactly. The associated knock-out (exchange) am- 
plitudes are the non-local effects in inelastic scattering. 
They contribute importantly to the evaluated scattering 
cross section, both in magnitude and shape pp. 

3 Results and Discussion 

3.1 MCAS and the low-energy regime 

In the paper detailing the MCAS method [5] cross sections 
and spectra from the scattering of neutrons from 12 C for 
energies to ~ 6 MeV were studied. The coupled-channel 
starting potentials were defined from a collective model 
which included quadrupole deformation. The spectrum of 
the target, 12 C, was truncated to just the lowest three 
states; the 0^ ground state, the 2+ state at 4.43 MeV and 
the 0jJ~ state at 7.96 MeV. We consider those results again 
in brief as they are the bound and scattering properties 
that the DPP we form will reproduce when used in a non- 
local Schrodinger equation. 

With parameter values as used previously [8j, MCAS 
calculations gave the energy variation cross section for low 
energy n+ 12 C elastic scattering that is shown by the solid 
curve in Fig. [TJ The filled circles are the data p!6p7] . The 
compound nucleus sub-threshold (bound) and resonance 
states that result are compared with the 13 C spectrum 
therein as well. The comparison of calculated results with 
data is very good at least to the energy (4.43 MeV) co- 
inciding with excitation of the first 2 + state. It is of note 
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Fig. 1. (Color online) MCAS results compared with data 
from the elastic n+ 12 C scattering and with the spectrum of 
13 C. Each 13 C state is identified by their parity and twice their 
spin. 



that, using MCAS for the negative (sub-threshold) energy 
regime gives the correct bound states of 13 C. Account- 
ing for the Pauli principle was crucial to find these re- 
sults; just how crucial was shown in recent publications [H 
[7]. Adding a Coulomb interaction also gave a good spec- 
trum for 13 N and good reproduction of fixed angle, proton 
elastic-scattering observables [5]. 

We have evaluated the DPP for this n+ 12 C system and 
results are shown in Fig. [2] Therein the Osi, Opi, and ds 

° 2 2 2 

wave DPP values are plotted in the top, middle and bot- 
tom segments for energies of 1.5 MeV and of 2.73 MeV 
in the left and middle panels respectively. The right panel 
contains some PB potential results that we discuss later. 
These DPP display a well shape with maximal depths on 
axis and in the nuclear surface region. The units for these 
DPP are MeV-fm -1 and the Osi wave contours have spac- 
ings of 20 MeV-fm -1 from a value of —20 MeV-fm -1 with 
the most outer (dashed) curves. The central well depth of 
these Osi wave DPP are -170 MeV- far 1 and -240 MeV- 

2 

fm -1 for 1.5 MeV and 2.73 MeV neutrons respectively. 
The 0p3 and Ods wave DPP respectively are given in con- 
tours starting with the outermost value of —5 MeV-fm -1 
with —5 MeV-fm -1 steps inward to maximal depth values 
of -22(-23) MeV-fm" 1 and -28(-32) MeV-fm" 1 for the 
incident energy of 1.5 (2.73) MeV. Clearly the DPP are 
strongly non-local and energy dependent; the 0s i wave for 
this case of n- 12 C scattering particularly so. The Opi and 
Ods wave forms are not as energy dependent over the to 
4 MeV projectile energy range but they, too, are markedly 
non-local. 

Values of DPPs along the diagonal (r = r') are shown 
in Fig. [3] The Osi and Op 3 wave DPP diagonal values, 
shown in the top and in the bottom (left) panels respec- 




Fig. 2. (Color online) Contours of the DPP found using 
MCAS for neutrons of energies 1.5 MeV (left) and 2.73 MeV 
(middle) compared with a PB energy-independent form (right). 




Fig. 3. (Color online) The radial variations of the DPP along 
the diagonal for Osi waves (top left), Op 3 waves (bottom left) 
and for Odj waves (right). Details are given in the text. 



tively, are for incident energies of 1.5 MeV (solid), 2.73 
MeV (long dashed), and 3.23 MeV (dashed curves). These 
clearly show the strong energy dependence of the 0s 1 wave 
interactions and the almost no energy variation in the 0p3 
wave DPP. The well shapes do not change with energy and 
the maximal depths are located at 2.0 and 2.2 fm respec- 
tively for the 0s 1 and Opa. wave results. In the right panel 
of Fig. El Odj wave DPPs for the energy of 2.73 MeV are 
displayed. The results for j = § and | are depicted by the 
dashed and solid curves respectively, illustrating the spin- 
orbit attributes formed using the MCAS approach. These 
variations were found [Sj without considering spin depen- 
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dent scattering data in the determination of the matrix of 
initial potentials that are input to the MCAS approach. 
Rather it was the spectrum (bound and resonant) of the 
compound nucleus ( 13 C) that set the parameterization. 
Nonetheless spin dependent scattering data were well pre- 
dicted by this model [SHE]. 



3.1.1 Results using the PB non-local potential 

As the non- locality of this model is energy independent, 
the only question is of what parameter values to use. As 
there arc no parameter values available that best fit rel- 
evant data, for simplicity we have chosen those used by 
Perey and Buck in their calculations with heavier nuclei. 
As noted, those parameter values suffice since as we seek 
a qualitative, but not quantitative, comparison of the fea- 
tures of the PB potentials with those found from MCAS. 
By using V NL = -74 MeV, R NL = 2.8 fm, a NL = 0.65 
fm and /3nl = 0.85 fm, calculations of Eq. ([20]) gave con- 
tour plots shown in the right panel of Fig. [2] The con- 
tour lines depict energy spacings of 10 MeV-fm -1 with 
the outermost contour that for —0.1 MeV-fm -1 . The cen- 
tral depths of these partial wave interactions are —46.1, 
-36.4, and -28.7 MeV-fm" 1 for the Osi, 0p 3 , and Ods 

2 2 2 

waves. Similarities exist between these non-local poten- 
tials and those found from low-energy MCAS calculations 
of the DPP. The depths of the wells for the two studies are 
comparable, with the depths of the DPP at 2.73 MeV be- 
ing close to the depths of the PB potential for the 0p| and 
Ods waves. The Osi wave MCAS depth is 3 times larger 

2 2 

than the PB value however. Nonetheless, it seems that the 
PB model does have non-locality qualitatively similar to 
what MCAS yields but the energy independence of the 
PB model's parameters seems too restrictive. 



3.1.2 Energy above the first excited state 

For energies where more than the elastic channel is open 
(above 4.43 MeV in the case of 12 C), there is the pos- 
sibility of flux loss to inelastic scattering. Therefore the 
DPP become complex. At 5.0 MeV in the n+ 12 C system, 
the MCAS method gives the complex DPP for the Osi 
and Opi waves that are displayed on the top and bot- 
tom in Fig. [1] respectively. Contours of the real and imagi- 
nary terms are shown on the left and right respectively in 
each case. The real parts of these potentials at this energy 
are repulsive with the central strengths of 5000 and 58.8 
MeV-fm _1 for the s and p wave cases respectively. The 
imaginary parts of both potentials are wells with minima 
of —620 and —8.6 MeV-fm -1 respectively. The contours 
for the real parts of these potentials are shown for every 
500 (10) MeV-fm- 1 with the 0s± (0p|) wave plots while 
the imaginary potentials for each indicate changes of 100 
(2) MeV-fm -1 . Clearly the forms of these potentials, as 
well as their being complex, have changed markedly from 
those found for nucleons with energies below that of the 
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Fig. 4. (Color online) Contour plots of the DPP for Osi 
wave (top) and Opi wave (bottom) at 5.0 MeV. The real and 
imaginary parts are shown on the left and right of each. 

threshold of the first excited state in 12 C. There is also 
distinctively new structure in the Opi DPP. 

The new structures of the DPP at 5.0 MeV are em- 
phasized with the plots of the diagonal values that are 
shown in Fig. [5] In the top panel, we show the Osi wave 
potential along the diagonal. The real part is depicted by 
the solid curve, the imaginary part by the dashed curve. 
The bottom section of the diagram contains the real and 
imaginary components of the Opi wave DPP along the 
diagonal. Again the real (imaginary) parts are given by 
the solid (dashed) curves therein. Note that the potential 
scales are shown in the brackets in each panel. Both DPP 
have large repulsive real parts with the Osi wave poten- 
tial being particularly strong. Both also have absorptive 
imaginary parts. The especial structure of the Opi case, 
already noted in Fig. [4] is most evident in this plot. 

Finally in Figs. [5] and [71 we show the energy variations 
of the diagonal (r = r') real parts of the DPPs. Those 
for the si-, di-, and d s -waves are depicted in the top, 
middle, and bottom of Fig. [5] respectively. Note that all 
of these strengths are given in units of MeV-fm -1 . Below 
the threshold (4.43 MeV) the potentials are purely real 
and attractive. There is quite strong energy variation in 
these DPP. At low energies, the DPP have the form of 
potential wells with minima in the nuclear surface region 
and those well depths increase with energy, most markedly 
as the energy approaches a value of about 4.8 MeV. For 
higher energies, those real potentials are strongly repul- 
sive with strength gradually decreasing. Passing through 
the threshold energy, the DPP acts like a potential wall 
at or about the nuclear surface. From the specification in 
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Fig. 5. (Color online) The diagonal values of the DPPs ob- 
tained using MCAS for 5.0 MeV neutrons scattering from 12 C. 
The Osi and Opi wave potentials are depicted in the top and 

2 2 

bottom sections respectively. Details are given in the text. 




r(fm) 

Fig. 6. (Color online) Energy variation of the DPP diagonal 
potentials for the si wave (top), for the d,3 wave (middle), and 
for the d& wave (bottom). 



Eqs. (fT5| and (TIT?]) , clearly the sharp change in the charac- 
ter of the real parts of these DPP do not necessarily occur 
at the threshold energy. Depending upon the sturmians 
and their eigenvalues the scale of change as well as the 
energy onset can vary with spin-parity of the scattering 
channel. With the latter, there seems to be a shift to an 



energy near 4.8 MeV for each case. However the size of 
change is very spin-parity dependent. 

The diagonal values of the real parts of the DPP for 
the pi and p| waves are depicted on top and bottom 
of Fig. [7] respectively. In both, the variation with energy 



V 200 




Fig. 7. (Color online) Energy variation of the DPP diagonal 
potentials for thepi wave (top) and for thepi wave (bottom). 



across threshold in the pi case is similar to that for the 

2 

d,3 wave though a double well aspect is evident now. The 
threshold effects upon the p3 wave DPP show changes but 
they are much less dramatic than in the DPP for the other 
partial waves. 



3.2 Non-locality effects for higher energy data 

The cross section and analyzing power for nucleons elasti- 
cally, and inelastically (to the 2 + , 4.43 MeV state), scat- 
tered from 12 C have been calculated using OBDME ob- 
tained from a complete (0 + 2)hui space shell-model calcu- 
lation [12]. Both WS and harmonic oscillator (HO) func- 
tions have been used with those (no-core shell-model) OB- 
DME to calculate the 12 C optical potentials. The tran- 
sition operators for both elastic and inelastic scattering 
processes are the medium dependent effective interactions 
at the relevant energies as described in Sect. 12.21 All of 
those elements are used in the DWBA98 program |15j . 
This distorted wave approximation (DWA) code includes 
exact evaluation of two-nucleon exchange amplitudes that 
define the knock-out process; amplitudes that are the re- 
sult of non-local aspects of the reactions since they involve 
the full nucleon density matrices of the target. 

With all details preset, just one calculation is made 
of both the elastic and inelastic observables for which 
the process then gives predictions. However, the elements 
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themselves need assessment of their quality. Much use [T] 
has established appropriate effective NN interactions in 
the nuclear medium as well as of the DWA method (when 
exchange amplitudes are treated exactly). Thus one need 
only assess, by other means if possible, the quality of the 
assumed structure. One excellent way to do that is to use 
the structure in analyses of electron scattering form fac- 
tors. 

3.2.1 An appropriate model of structure 

For the specific case we study, 12 C, the longitudinal form 
factors for elastic and inelastic (0 + — * 2 + ) reactions, and 
the transverse electric form factor with the latter, have 
been measured accurately over a reasonable range of mo- 
mentum transfer values. Those form factors are compared 
with our calculated ones in Fig. [51 Considering the elastic- 




q (fm" 1 ) 

Fig. 8. (Color online) Electron scattering form factors for 12 C. 
The elastic longitudinal form factor is shown in the left panel 
while the inelastic ones (from excitation of the 2 + (4.43 MeV) 
state are shown on the right; longitudinal (top) and transverse 
electric (bottom). All details are given in the text. 

scattering results (left panel) first, the data are those of 
Jansen et al. 2D] (squares), of Sick and McCarthy |21 (di- 
amonds), and of Nakada et al. [2% (circles). The dashed 
and solid curves are the results found using WS single 
particle wave functions with shell occupancies from OTuo 
and (0 + 2)fiuj shell-model calculations respectively. The 
WS potential parameters used to specify the bound-state 
wave functions and the shell occupancies are given in the 
review pQ. Similar results were found using HO single- 
nucleon wave functions when the oscillator length was 1.6 
fm. All results agree well with the data though the higher 
momentum transfer values favor the larger space model of 
structure. But the form factors for the excitation of the 
first excited (2 + ) state given in the right panel of Fig. [8J 
differentiates more strongly. 

In the top part of this panel, we present data and re- 
sults for the longitudinal electric form factor. The data 



displayed by the open squares and filled circles are two sets 
reported by Flanz et al. [23j . They are compared with re- 
sults found using shell-model transition OBDME from the 
Ohuj (dashed curve) and from the (0 + 2)hu (solid curve) 
shell-model calculations. Clearly the additional contribu- 
tions from transitions out of the Op shell that result in the 
larger space shell-model study give the extra transition 
strength required to match observation. 

The solid curve shown in the bottom panel is the trans- 
verse electric form factor calculated using the (0 + 2)huo 
shell-model transition OBDME and WS wave functions. 
The match to data [23] is good especially when one notes 
that the separate proton and neutron contributions to 
this form factor (shown by the dash-double dotted and 
dot-double dashed curves respectively) have amplitudes 
that interfere destructively to determine the total result. 
In the definition of the transition form factor, Siegert op- 
erators were used [I] to account for meson exchange cur- 
rent corrections. Thus we have confidence that the no-core 
(0 + 2)Twj shell model describes well the ground state of 
12 C and the excitation of its first excited state, as well as 
of using those wave functions and OBDME in analyses of 
nucleon scattering from 12 C. We consider two cases; those 
for 95 MeV neutrons and for 200 MeV protons. 

3.2.2 Credibility of the effective interaction 

The differential cross sections and analyzing powers for 
elastic scattering of 200 MeV protons from 12 C that result 
on using the full non-local optical potentials are compared 
with data [24] in Fig. [9] Consistent with the results found 
with the elastic electron scattering form factors, using the 
Cohen and Kurath [55] Op shell-model structure makes 
but little change to the proton elastic-scattering calcula- 
tions from those displayed. In the left and middle panels 
in this figure, the solid and long dashed curves display the 
results obtained by using WS and HO bound-state wave 
functions. Those panels contain results found from calcu- 
lations made using t- (left) and g- (middle) folding. In t- 
folding, the purely free two-nucleon i-matrices were used 
as the effective interaction. Thus comparing the results 
of these two panels shows how important are the medium 
modifications of the effective NN interaction. It is evident 
that ^-folding with WS bound states gives an optical po- 
tential from which a differential cross section is obtained 
that matches the data best. The importance of medium 
modification in the effective NN interaction is even more 
obvious with the analyzing power. The <?-folding results 
are in much better agreement with data than are the t- 
folding ones. However, with this observable as well, there 
is little to choose between the results obtained using HO 
and WS bound states for the single-particle bound states 
of 12 C. 

The effect of omitting the exchange amplitudes in defin- 
ing the elastic scattering are shown in the rightmost panel. 
Therein the solid curves are the complete (/-folding model 
results also shown in the middle panel, while the dot- 
dashed curves are the results when the exchange ampli- 
tudes are ignored. Clearly those amplitudes are essential 
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Fig. 9. (Color online) 200 MeV proton elastic-scattering cross 
section (top) and analyzing power (bottom). Data [24\ (circles) 
are compared with the results found using the t- (left) and g- 
folding (middle) optical potentials The right panel shows the 
effect of omitting the knock-out amplitudes of the gr-folding 
potentials. All details are given in the text. 



in finding good predictions for both cross sections and 
analyzing powers. Further discussion is given in the next 
sub-section. 




of 12 C. Data for this reaction are compared with micro- 
scopic DWA model results in Fig. [TUJ Comparison of the 
two results with the data demonstrates the important role 
played by medium effects in the AW interaction in this 
process. The variation in shapes of results for scattering 
angles above 30° in the center of mass indicates that. In 
addition the strength of the scattering is well predicted 
by the (0 + 2)7iu; structure model. When a simple OTilu 
spectroscopy was used, this cross section was a factor of 
four too weak in comparison with the data [I] requiring 
upward scaling equivalent to an effective charge of 0.5e to 
match measurement. 



20 40 60 80 3.2.3 Effects of the non-local exchange amplitudes 



As the non-local nature of the optical potentials formed 
by g-folding can be inconvenient, in the past the exchange 
terms from which those non-localities arise either have 
been ignored or approximated by a local potential. While 
the latter approach is the more sensible, we consider what 
happens by ignoring the exchange terms and so retaining 
only the "direct" potentials; the leading terms of Eq. (|2"o| 
i.e. those formed by folding with just the nuclear density 
and not the full density matrices. By comparing complete 
with pure direct calculation results, the influence of the 
exchange terms can be defined. 

In the right-most panel of Fig. [9j the differential cross 
section and analyzing power data (200 MeV elastic proton 
scattering) are compared with the results of calculations 
made with the exchange amplitudes included, without any 
simplification, and with those exchange amplitudes totally 
excluded. The results are shown by the solid and dot- 
dashed curves respectively. The difference between the two 
sets of results is as noteworthy as the quality of agreement 
between the full calculation results and the data. Such dis- 
parity between results with and without the exchange am- 
plitudes persist over a wide range of energies pQ . Of par- 
ticular note is that the direct and exchange contributions 
destructively interfere to produce the final result. Also the 
two amplitudes lead to very different momentum-transfer 
effects in order that the cross-section results end up in the 
agreement with data found. 
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Fig. 10. The cross sections and analyzing powers for the 
inelastic scattering of 200 MeV protons exciting the 2 + 4.43 
MeV state of 12 C. The solid curves are the complete results and 
the dashed curves are those found when the medium effects in 
the effective NN interaction are omitted. 



The quality of the structure and the import of us- 
ing medium modified effective NN interactions are con- 
firmed by the results obtained for the inelastic scattering 
of 200 MeV protons exciting the the 2+ (4.43 MeV) state 



3.2.4 Neutron scattering and effects of exchange amplitudes 

The cross section from the elastic scattering of 95 MeV 
neutrons from 12 C has been measured recently |26j and 
these data are shown by the filled circles in the left panel 
of Fig. HU Older data taken at 40.3 MeV [27] are shown 
therein by the filled squares. The middle and left pan- 
els show the results of 95 MeV neutron scattering elastic 
(top) and inelastic, to the 2 + (4.43 MeV) state, (bottom) 
cross section and analyzing powers respectively. The solid 
curves depict the results obtained using the complete in- 
folding model with WS wave functions and orbit occupan- 
cies from the (0 + 2)hu> shell model. The dashed curves 
depict what result on omitting the exchange amplitudes. 
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Fig. 11. In the left panel are the cross sections for 40.3 MeV 
and for 95 MeV neutrons scattered elastically from 12 C. In 
the middle panel are the cross sections for elastic (top) and 
inelastic scattering (bottom) found with (solid) and without 
(dashed) the knock-out amplitudes. The right panel are the 
analyzing powers associated with them. Details are given in 
the text. 

At both energies, treating the exchange terms with- 
out approximation to get the optical potentials gave cross 
sections in very good agreement with the data. The con- 
trast of these results with those given by the calculations 
made omitting the exchange amplitudes is dramatic. The 
structure of the cross sections and analyzing powers are 
radically changed as is the predicted magnitude of the 
inelastic excitation. These variations are similar to those 
noted above with the (elastic) scattering of 200 MeV pro- 
tons. 

Of particular import is that, at all three energies con- 
sidered, the role of the exchange (non-locality) in eval- 
uations of both the elastic and inelastic-scattering cross 
sections is of destructive interference. Such interference 
makes representation of the non-local interactions by an 
equivalent local interaction suspect. Also the effects in the 
analyzing powers are dramatic and it must be remembered 
that those observables involve the cross section as a de- 
nominator. A fit to any analyzing power then without first 
having a good cross section prediction is but fortuitous. 



4 Conclusions 

The origins and effects of non-localities in the NA opti- 
cal potential have been studied with separate methods of 
data analysis most appropriate for low and for medium 
projectile energies. The scattering of nucleons from 12 C 
have been considered. 

For low energies (0 to ~ 6 MeV) the MCAS method 
was used to describe n+ 12 C scattering. With that theory 
the origins of non-localities in optical potentials due to 
antisymmetry of projectiles with target nucleons and the 
dynamic polarization potential representation of channel 
coupling were detailed. A collective model was used to 
specify the matrix of interaction potentials that were the 



basic input to the approach. The resulting DPP is strongly 
non-local with also strong /-dependence. At the energies 
studied the dominant term is that of 0s i wave interaction 

2 

whose character varies markedly to carry the resonance 
features of the full coupled-channel results. 

At medium energies, the g-folding method with no- 
core, large space, Shell model wave functions and the Mel- 
bourne force, a complex effective NN interaction that is 
density and energy dependent, was used. That approach 
attributes the effective channel coupling to be with an 
effective infinity of target states in which one or more nu- 
cleons are in the continuum. It is equivalent to using the 
KMT theory but with the effective interaction between the 
projectile and each and every target nucleon being mod- 
ified by the medium. The cross sections obtained when 
compared with data, established that one cannot ignore 
either medium modification of the NN interaction or the 
exchange (knock-out) contributions in forming optical po- 
tentials. That is observed strongly also when the DWA 
is used to evaluate observables from inelastic scattering. 
The complete <?-folding process makes the optical poten- 
tial complex, energy dependent, and non-local. But, as 
the ^-matrices are strongly medium dependent, the non- 
local attributes of the optical potentials themselves are 
not easily displayed. Indeed in the DWBA98 code such 
are not explicitly evaluated in finding solutions of the 
integro-differential equations from which phase shifts are 
specified pQ. The import of that non-locality is evident 
however in the comparisons made with data at many en- 
ergies and with and without the exchange terms included 
in the calculations of the cross sections and analyzing pow- 
ers. Clearly, any localization of the non-local potential is 
approximating a most important factor in data analysis. 
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